Generation of Large Moments in a Spin-1 Chain with Random Antiferromagnetic 
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We study the spin-1 chain with nearest neighbor coupUngs that are rotationally invariant, but 
include both Heisenberg and biquadratic exchange, with random strengths. We demonstrate, using 
perturbative renormalization group methods as well as exact diagonalization of clusters, that the sys- 
tem generates ferromagnetic couplings under certain circumstances even when all the bare couplings 
are antiferromagnetic. This disorder induced instability leads to formation of large magnetic mo- 
ments at low temperatures, and is a purely quantum mechanical effect that does not have a classical 
counterpart. The physical origin of this instability, as well as its consequences, are discussed. 
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Classical spin models have been studied extensively 
over the past half century. The inclusion of quantum 
mechanical nature of the spin variables result usually in 
quantitative effects such as renormalization of transition 
temperatures, order parameter in the ordered phase, spin 
wave velocity, etc. However, it began to be realized about 
two decades back that quantum fluctuations could result 
in qualitative changes such as the nature of the phase it- 
self. Ma et al. [|l[ showed that in a spin- 1/2 chain with 
nearest neighbor antiferromagnetic Heisenberg interac- 
tions of random strength, quantum fluctuations lead to a 
divergent density of excitations at low energy scales even 
in the absence of such divergence in the distribution of 
bare couplings. Around the same time, Bhatt and Lee 
1^ independently showed that this occurs for highly dis- 
ordered spin-1/2 antiferromagnets with short range in- 
teractions in higher dimensions as well, and that this 
phenomenon is a purely quantum mechanical effect [pj. 
A few years later, Haldane showed that the uniform 
integer spin chains with antiferromagnetic Heisenberg in- 
teractions exhibited a gap, and concomitant short-range 
spin-spin correlations in the ground state, in contrast 
with the half integer as well as classical spin chains. The 
Haldane gap, which scales as exp(— ttS') with the spin S, 
is also a purely quantum mechanical effect, unobtainable 
from the classical {S — > oo) limit, where it appears as 
an essential singularity. In the past few years, studies 
of the quantum Ising model in a transverse field in one, 
two and three dimensions [gJI] as well as Heisenberg spin 
chains RM with randomness, have demonstrated that 
the ground state (T = 0) phase diagram includes a Grif- 
fiths phase with divergent response functions due to rare 
fluctuations. Thus, electronic systems in the quantum 
regime display a richer variety of phenomena resulting 
from the interplay between correlation and disorder, than 
is usual for classical systems. 

In this paper, we describe another phenomenon occur- 
ring in spin-chains as a purely quantum mechanical effect 
- namely, the generation of ferromagnetic (F) couplings. 



and consequently large moments leading to a Curie sus- 
ceptibility at low temperatures, in a spin-1 chain with 
isotropic but random antiferromagnetic (AF) couplings. 
Though the generated moments become arbitrarily large 
in the low temperature limit, the phenomenon of the 
generation of ferromagnetic couplings relies on a purely 
quantum mechanical effect, and has no classical analog. 

For the random spin-1/2 chain with near neighbor AF 
interactions, Fisher [Q showed that the real space renor- 
malization group (RG) scheme |Qj|] becomes asymptot- 
ically exact, and leads to a "random singlet" phase, 
where distant pairs of spins form singlets in a hierar- 
chical manner dependent on the realization of the ran- 
dom bonds, and dominate the low energy physics. In the 
presence of randomly placed ferromagnetic couplings of 
arbitrary concentration, however, Westerberg et al. pd| ] 
showed that the random singlet phase is destroyed in 
one-dimension due to the formation of large moments by 
active, ferromagnetically coupled spins, and the magnetic 
susceptibility at asymptotically low temperature assumes 
a pure Curie (1/T) form, right upto the ferromagnetic 
point at zero concentration of antiferromagnetic bonds. 

For a spin-1/2 system, spin rotational symmetry 
uniquely constrains the coupling between spins to be 
the Heisenberg form ||l^, JSi ■ Sj. For spins with 
S > 1/2, however, the most general form of isotropic cou- 
pling between spins i and j with spin rotational invari- 
ance may be written in powers of the Heisenberg form: 
Hij — J2n=i J^^^H^i ' Sj)"- Foi' the random spin-1 chain 
with nearest neighbor couplings, this implies a Hamilto- 
nian written most generally as: 



H — ^^[JiSi ■ S,;+i + Di{Si 



S.+i)'] 



^Y.yJ^+ A'[cose.S, • S,+i -f sin0,(S, . S,+i)2], (1) 

i 

where the J's and Z3's are uncorrelated random variables. 

With purely Heisenberg couplings, it was showed H] 

that as long as there are no F bonds in the bare Hamil- 



tonian, the system cannot be in the pure Curie paramag- 
netic phase. In this paper we show that this is no longer 
the case when biquadratic couphngs are present u%. 

With no randomness (i.e. Ji — J,Di = D), the prop- 
erties of the Hamihonian of Eq. (2) is controUed by the 
angular variable 6 satisfying tan9 = D/J. There exist 
four different phases (see Fig. |l|). For 7r/2 < 6 < 57r/4, 
each individual bond favors a total spin Stot — 2 state for 
the pair it connects, and the ground state of the entire 
chain is the spin fully polarized ferromagnetic state. For 
— 7r/4 < 9 < 7r/4 the system is in the Haldane gapped 
phase B. For — 37r/4 < 6 < — 7r/4 the chain is sponta- 
neously dimerized |14|, while an extended gapless phase 
has been predicted in the region 7r/4 < < 7r/2 [[l5[ . 
Except for the ferromagnetic phase, the other phases all 
have singlet ground states with unbroken spin-rotational- 
symmetry, and we therefore refer to bonds in this region 
as antiferromagnetic. 

In the presence of strong randomness, the system may 
be studied using a hierarchical real space RG approach. 
We search for the bond in the system with the largest 
gap separating its ground state and lowest energy excited 
state, say the bond coupling spins 2 and 3, with coupling 
constants J2 and D2 (see Fig. 2a). If the ground state 
of this bond is a singlet (— 37r/4 < 62 < arctan^), then 
spins 2 and 3 form an inert singlet in the low-energy 
states of the system, and mediate effective couplings be- 
tween their neighboring spins 1 and 4, which may be 
calculated using second order perturbation theory Mn : 



Ji. 



D, 



(2Ji-gi)(2J3-A3 

3( J2 - 3-D2) 
2D1D3 



D1D3 



HJ2-D2)' 



9( J2 - D2) 



(2) 
(3) 



If the ground state of the bond is a triplet (arctan i < 
02 < 7i'/2), spins 2 and 3 form an effective spin 2' with 
5*2' = 1, and its couplings to its neighbor spin 1 are 



Jl2' 
Dl2' 



(Ji-i?i)/2; 
-Di/2. 



(4) 
(5) 



Couplings to spin 4 have identical expressions. Should 
the ground state of the bond be a quintuplet, an effective 
spin with 5 = 2 forms, and the structure of the original 
spin-1 chain gets distorted. 

Examination of Eqs. (2) and (0) shows that the gen- 
erated bond may be ferromagnetic, even if all the bonds 
involved are antiferromagnetic and favor singlet ground 
states. For concreteness, we consider the case where 
bonds 2 and 3 are Heisenberg, i.e., D2 = -D3 = 0. In 
this case the effective bond between spins 1 and 4 are: 
Ji4 = 2J3(2Ji - Di)/3J2 and D^ = 0. Therefore the 
generated bond is Heisenberg and ferromagnetic if 



Ji - Di/2 < 0. 



(6) 



Clearly, bond 1 can be AF and satisfy Eq. (|^) if it lies 
in the shaded region of Fig. ^ 

We have verified the above results of perturbation the- 
ory, that it is possible to get Stot = 2 ground states when 
Eq. (H) is satisfied, by performing exact diagonalization 
of four spin clusters. The results of one such calculation 
is shown in Fig. Bt as a function of J2, the ground state 
of the cluster changes from Stot = (which evolves from 
a product of singlets between spins Si & S2 and S3 & ^4 
in the limit J2 = 0) to Stot = 2. 

The perturbative RG described above is reliable when 
the randomness is strong and the distributions of the cou- 
plings are broad [|l],|0|. We now consider the opposite 
limit of dilute randomness, namely a uniform AF spin-1 
chain in the Haldane phase ("7r/4 < 6 < 7r/4, see Fig. 
|l|), with a finite Haldane gap A and a small fraction of 
impurity bonds that are much weaker than A (see Fig. 
2b). Here one must identify the true low energy degrees 
of freedom ||^. If the impurity bonds were taken away, 
the original chain would have been chopped into decou- 
pled segments; the low-energy degrees of freedom are the 
two spin- 1/2 at the two edges of each segment with a 
coupling (that can be either F or AF) decreasing expo- 
nentially with the length of the segment. Putting back 
the impurity bonds, as long as they are weak compared 
to A, does not alter the bulk structure of the segments; 
their primary effect is to couple neighboring edge spins in 
different segments. Let us assume bond 1 coupling spins 
1 and 2 is such an impurity bond, with Ji,Di <^ A. To 
calculate the coupling between the two edge spin-l/2s, 
which we label 1' and 2', we project the original opera- 
tors onto the subspace of states below the Haldane gap, 
i.e., states of the effective edge spins [ p7[ : 

Hv2' = PH12P = JiPSi • S2P + i:'iP(Si • S2)2p, (7) 

where P is the projection operator. Rotational sym- 
metry as well as properties of spin-1/2 guarantee that 
Hi'2' = Ji'2'Si' • S2' + Ci'2', where Ci'2' is a constant. 
We also have PSi • S2P = PSiP • PS2P, because spins 
1 and 2 live in decoupled Hilbert spaces if bond 1 were 
not there, and the Wigner-Eckart theorem jl^ guaran- 
tees PSjP — aS^. The constant a depends on bulk 
properties of the segments; for infinitely long segment 
with Heisenberg coupling, a w 1.0640 |1^. Similarly 
P(Si ■ S2)'^P — a + 6S1' • S2', and the important con- 
stant b may be determined by calculating certain ma- 
trix elements of (Si • 82)^ in the subspace: b = 2[('fi/|2' 

|(Sl • S2)2| Tl't2') - (Tl'i2' |(Sl • S2)2| tl'i2')]. Us- 

ing the commutation relations of S and the fact M] 
(Ti' \Si\ Ti') = a/2 etc, we obtain 6 = ~a^/2. There- 
fore Ji'2' = a^{Ji — Di/2), and again we find Ji'2' rnay 
be ferromagnetic even when the original bond 1 was AF, 
with the same condition as before, namely Eq. (o). We 
have thus demonstrated the existence of such an insta- 
bility to form ferromagnetic couplings in both the high 
disorder and dilute disorder limits. 



We now discuss the origin of this ferromagnetic insta- 
bihty, and in particular, the significance of the special 
combination of J and D in Eq. (g). By introducing a 
new coupling constant K — J — -D/2, the RG equations 
(2) and (4) simplify significantly: 



K^ 



4K1K3 



3(^2 - 5D2/2) 



K, 



Ki/2. 



(8) 



Combining these with equations (||) and (^), we find K 
and D decouple except through energy denominators, and 
one cannot generate K from D, or vice versa. This sug- 
gests that K and D represent couplings of operators with 
different symmetry properties. 

To proceed further, we note that products of differ- 
ent components of the spin operator (that appear in the 
Hamiltonian) may be organized to form traceless irre- 
ducible spherical tensor operators ||l^: Yim{S), which 
is defined by replacing cos 9 by Sz, sin 9 cos </> by Sx 
and sin 9 sin by Sy in the usual spherical harmonics 
Yim{9, 0), and symmetrizing noncommuting components. 
^im(S) = for I > 2S. A general way to write down the 
coupling between two spins that respect rotational sym- 
metry (equivalent to Hij = J2n=i J^^'H^i ' ^j) 



IS 



2S 



i/y^^xC) ^ (-f)™il,„(sorz,_,„(s,), (9) 



1=0 



=-i 



where K^ is the coupling constant of rank I spherical 
tensors. It is easy to verify that 

. 1 

Si • S2 = y Y. (-l)'">lm(Sl)yi,_,„(S2), (10) 

rn— — 1 



27r 



(Si • 82)2 = -4^ ^ (-i)™ri„(Si)yi,_™(S2 



m — — 1 



8 2 

+ if E (-i)"^2™(Si)y2,-™(S2) + -s?s2. (11) 

m=-2 

Therefore for the coupling of the form JSi • S2 + -D(Si ■ 
S2)^ we have iffi' = ^i2J-D) oc if , and K^^) = ^D. 
Thus the K variable is proportional to the coupling be- 
tween rank 1 tensors (vectors), nothing but the Heisen- 
berg coupling (in this tensor representation) ; and Eq. m) 
indicates the Heisenberg coupling is ferromagnetic. 

The advantage of writing the Hamiltonian in terms 
of couplings of irreducible spherical tensors instead of 
powers of Si ■ S2 is that different symmetry properties of 
tensors with different ranks do not allow them to mix in 
1st and 2nd order perturbation calculations, as we have 
already seen; while (Si • S2)" in general includes tensor 
couplings with all ranks up to n. For example, we know 
yim(S) acting on a singlet creates an eigenstate of Sj^j 
and S^g^ with Stot = I and S^^^ = m, therefore in the 
second order perturbation a coupling of rank / between 



spins 1 and 4 is mediate through the channel of excited 
states with Stot = I of spins 2 and 3. More remarkably, 
when projecting spin-1 couplings to couplings between 
spin-1/2 edge spins, the fact that spin-1/2 object does 
not support tensors with ranks higher than 1 guarantees 
that the coupling must be proportional to J — D/2 of the 
original coupling, and the original rank 2 tensor coupling 
simply gets eliminated. 

The above discussion also gives us insights into the 
origin of the ferromagnetic instability we demonstrated. 
Even though a single bond favors a singlet ground state, 
it may well contain ferromagnetic couplings between ten- 
sors with certain rank (say, rank 1, or Heisenberg cou- 
pling). In the absence of randomness no ferromagnetic 
instability is triggered by such couplings. In the presence 
of randomness, however, energy scales in the system gets 
separated, and low- and high- energy subspaces of the 
Hilbert space get perturbatively decoupled (which is the 
basis for perturbative RG). In projecting to low-energy 
subspaces, certain AF couplings may get suppressed for 
symmetry reasons (as we have seen), while the original 
subdominant ferromagnetic couplings may survive and 
become dominant. In the example we illustrated above, 
the AF rank 2 couplings are suppressed in the low energy 
subspace, while the ferromagnetic rank 1 (Heisenberg) 
coupling survives. This is the origin of the ferromagnetic 
instability, which is a purely quantum mechanical effect, 
and absent without randomness. 

With the generation of ferromagnetic bonds between 
segments as well as within segments, the decimated chain 
is very different from that obtained with pure Heisenberg 
coupling by Hyman and Yang M. In fact, it becomes of 
the universality class of the spin-1/2 chain with random 
AF and F couplings studied by Westerberg et al. [O. 
Using our formalism, their RG scheme may be easily 
generalized to include higher spins and non-Heisenberg 
couplings, as flows of couplings of different rank tensors 
tend to decouple. Our results show that at low ener- 
gies high rank couplings are strongly suppressed com- 
pared to Heisenberg couplings, if the original couplings 
are dominantly Heisenberg, which is a likely situation 
in nature |20|. Therefore the active degrees of freedom 
in our case will be weakly coupled large moments as in 
the case of Westerberg ei al., which lead to a pure Curie 
susceptibility in the low temperature limit, because of a 
perfect cancellation between growing moments /^ and re- 
ducing numbers N in the expression for the susceptibility 
X = N{T)n{TflT as T is lowered. 

In summary, we have demonstrated the possibility 
of generating ferromagnetic couplings in quantum spin 
chains with random rotationally invariant couplings, 
starting from the antifcrromagnctic sector of the phase 
diagram of the pure chain. Though we have discussed 
S — 1 exclusively for concreteness, similar effects would 
be expected for higher spins, though their phase diagrams 
even for the pure case are not fully known. For integer 



spins, our arguments are valid in both large disorder and 
dilute limits, and the former case applies to larger half- 
odd-integer spins as well. However, as our treatment 
shows, the effects are purely quantum mechanical, and 
dependent on either the spin gap or the region of valid- 
ity of second order perturbation theory, both of which 
are known to become smaller with increasing S. Conse- 
quently, we expect the region of this anomalous behavior 
to decrease with increasing S, and disappear in the clas- 
sical limit S* — *■ oo. 
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FIG. 1. Phase diagram of a pure spin-1 chain. Solid lines 
are phase boundaries. Shaded regions represent couplings in 
the antiferromagnetic sector satisfying J — D/2 < 0, which, 
in the presence of randomness, could generate effective ferro- 
magnetic bonds at low energies. 



FIG. 2. Illustration of spin decimation procedures, (a) 
Strong randomness case. When the strongest bond (between 
spins 2 and 3) has a singlet ground state, spins 2 and 3 are 
decimated and an effective bond connecting 1 and 4 is gen- 
erated; when the ground state is a triplet, 2 and 3 form an 
effective spin-1 object ^2', which is coupled to its neighbors 
1 and 4. (b) Dilute randomness case. When a uniform spin-1 
chain in the Haldane phase is broken into finite segments cou- 
pled by weak impurity bonds (dotted lines), the low energy 
degrees of freedom are half spins living at the edges of each 
segments; there is weak coupling between neighboring half 
spins, both in the same segment (broken lines) and different 
segments (dotted lines). 
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FIG. 3. Ground state level crossing in a four-spin cluster 
for Ji = -1.0, Di = -1.5, Ja = 1.0, D2 = D3 = 0. The plot 
shows energies of the lowest energy quintuplet {Stot = 2, solid 
line) and triplet {Stot = 1, broken line) states, measured from 
lowest energy singlet state. For J2 > 3.46, the ground state 
has total spin Stot = 2, despite the fact that all three bonds 
favor singlet ground states. 
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